Atomic CP-violating polarizability 
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Searches for CP violating effects in atoms and molecules provide important constrains on com- 
peting extensions to the standard model of elementary particles. In particular, CP violation in an 
atom leads to the CP-odd (T,P-odd) polarizability f3 C¥ ': a magnetic moment /i CP is induced by an 
electric field So applied to an atom, /i CP = f3 CF £o- We estimate the CP-violating polarizability for 
rare-gas (diamagnetic) atoms He through Rn. We relate f3 CF to the permanent electric dipole mo- 
ment (EDM) of the electron and to the scalar constant of the CP-odd electron-nucleus interaction. 
The analysis is carried out using the third-order perturbation theory and the Dirac-Hartree-Fock 
formalism. We find that, as a function of nuclear charge Z, /3 CP scales steeply as Z 5 R(Z), where 
slowly-varying R(Z) is a relativistic enhancement factor. Finally, we evaluate a feasibility of set- 
ting a limit on electron EDM by measuring CP-violating magnetization of liquid Xe. We find that 
such an experiment could provide competitive bounds on electron EDM only if the present level of 
experimental sensitivity to ultra- weak magnetic fields [Kominis et al, Nature 422, 596 (2003)] is 
improved by several orders of magnitude. 

PACS numbers: 11.30.Er,32.10.Dk,31.30.Jv 



I. INTRODUCTION 

The existence of a permanent electric dipole moment 
(EDM) of a particle simultaneously violates two discrete 
symmetries: parity (P) and time reversal (T). By the 
virtue of the CPT-theorem, the T-violation would imply 
CP-violation P, 0. While no EDMs have been found 
so far, most supersymmetric extensions of the Standard 
Model of elementary particles predict electron EDMs, d e , 
that are within a reach of planned and on-going exper- 
imental searches. Here we investigate a related T-odd, 
P-odd quantity — CP- violating polarizability, (3 CP , in- 
troduced recently by Baryshevsky 3] . For a diamagnetic 
atom, a non-vanishing (3 CP could provide an unambigu- 
ous signature of the electron EDM or other CP-violating 
mechanisms. Here we relate /3 CP to d e via ab initio rela- 
tivistic calculations for closed-shell atoms. We also relate 
(3 CP to the scalar constant of the CP-odd electron-nucleus 
interaction. 

An interaction of an atom with external DC electric 
field in the presence of the electron EDM causes spin 
polarization in the direction of the field |4j . The first at- 
tempt to measure corresponding magnetization of the fer- 
romagnetic crystal was made by Vasiliev and Kolycheva 
in 1978 5]. According to Lamoreaux 6], modern tech- 
niques allow to improve that old measurement by many 
orders of magnitude and reach the sensitivity, which al- 
lows to improve present limit on the electron EDM 0: 

d e (Tl) < 1.6 x lCT 27 e-cm. 

First results of the new generation of experiments with 



ferromagnetic solids were recently reported by Hunter j^] • 
Characteristic feature of the experiments with macro- 
scopic magnetization is the dependence of the signal on 
the density of atoms. That gives a huge enhancement in 
sensitivity for a condensed phase sample. 

R is generally assumed that diamagnetic atoms are not 
useful for the search of the electron EDM. However, Bary- 
shevsky has recently pointed out Q that CP-violating 
magnetization would also exists for a diamagnetic atom. 
For a spherically-symmetric atom, the E-field-induced 
magnetic moment /z CP can be expressed in terms of CP- 
violating polarizability /3 CP as 



CP 



= {3 CP £o 



(1) 
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where £o is the strength of the electric field. This ob- 
servation opens new experimental possibilities. For ex- 
ample, one can measure magnetization of liquid xenon 
in a strong external electric field. The advantage of the 
experiment with diamagnetic liquid in comparison to fer- 
romagnetic solids is a much lower magnetic noise. 

For a diamagnetic (closed-shell) atom the magnetiza- 
tion appears in the higher orders of the perturba- 
tion theory than for the open-shell atoms. In this paper 
we calculate polarizability j3 for rare-gas atoms He 
through Rn using third-order perturbation theory and 
Dirac-Hartree-Fock (DHF) formalism. 

Further, we evaluate a feasibility of setting a limit on 
electron EDM by measuring CP-violating magnetization 
of liquid Xe (LXe). We consider the effect of the envi- 
ronment on /3 CP of Xe atoms in LXe. We use a simple 
cell model of an atom confined in a spherically-symmetric 
cavity 9j. In a non-polar liquid, such a cavity roughly 
approximates an averaged interaction with the neighbor- 
ing atoms. We solve the DHF equations with proper 
boundary conditions at the cavity radius. For LXe, we 
find that compared to the CP-odd polarizability of an 
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isolated atom, the resulting CP-odd polarizability of an 
atom of LXe is suppressed by about 65%. 

We find that the CP-violating polarizability exhibits 
an unusually strong de pen dence on the nuclear charge 
Z . Previously, Sandars [lflllll ] has shown that an atomic 
enhancement factor for the electron EDM is of the order 
of a 2 Z 3 , where a « 1/137 is the fine structure constant. 
As we demonstrate below, for a diamagnetic atom, the 
polarizability /3 CP vanishes in the non-relativistic approx- 
imation. Because of that it is suppressed by a factor of 
(aZ) 2 . With the Sandars' enhancement factor this leads 
to a steep, Z 5 , scaling of the effect. 

Recently there was a renewed interest to the CP-odd 
weak neutral current interactions of electrons with nu- 
cleons |T^ . It is known that in atomic experiments the 
electron EDM is indistinguishable from the scalar CP- 
odd weak neutral currents Q, Il3| . Any new limit on 
the electron EDM from the atomic experiments will also 
lead to the improved limit on the scalar constant of the 
CP-odd electron-nuclear interaction. Here we relate com- 
puted (3 to the scalar constant of the CP-odd electron- 
nucleus interaction. 

The paper is organized as follows: In Section [H] we de- 
rive the third-order expression for the CP-violating po- 
larizability and use the independent-particle approxima- 
tion to simplify the atomic many-body expressions. In 
Section IIIII we present results of our DHF calculations 
of (i CP for rare-gas atoms and derive the Z-scaling of 
/3 CP . In Section El we evaluate a feasibility of setting a 
limit on electron EDM by measuring CP-violating mag- 
netization of liquid Xe. Finally, in Section [V] we draw 
conclusions. Unless specified otherwise, atomic units 
|e| = h — m e = 1 and Gaussian system for electro- 
magnetic equations are used throughout. In these units, 
the Bohr magneton is fis = a/2 and the unit of magnetic 
field is m^/ft 4 « 1.72 x 10 7 Gauss. 



II. FORMALISM 



tion III CI For a spherically-symmetric system, the in- 
duced magnetic moment will be directed along the ap- 
plied E-field. 



A. Third-order formula for the induced magnetic 
moment 

We develop the perturbative expansion for the atomic 
wavefunction |^o) m terms of the combined interac- 
tion W — V CP + T/ cxt . Here V cxt is the interaction 
with the external electric field applied along the z-axis, 
yext _ _£) z £ 0j jy z being the z-component of the elec- 
tric dipole moment operator. To estimate the dominant 
contribution to (fj,), it is sufficient to truncate the per- 
turbative expansion for the atomic wavefunction at the 



second order in W, l^o) 



I^ 0) ) + K' 



I* 



(2)\ 



Then 



the expectation value of the magnetic moment reads 

To arrive at the above expression we used a simplify- 
ing fact that the magnetic moment is a P-even operator, 
while both | \I/q ^ } and I^q 2 "*) have parities opposite to the 
one of the first-order correction I^q 1 ^). 

The textbook expressions for the first and second-order 
corrections to wavefunctions can be found, for example, 
m Ref. 0. With these expressions, 



<M CP > - <M CP ) 1 + <M OP > 9 + <M CP ), > (3) 

(4) 
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(5) 



(6) 



In this Section we derive the expression for CP- 
violating polarizability within the third-order perturba- 
tion theory. Further, we simplify the derived expression 
using the Dirac-Hartree-Fock approximation for atomic 
many-body states. 

The problem to be solved can be formulated as fol- 
lows: What is the induced magnetic moment (fi) of an 
atom perturbed by an external electric field £$1 It is 
easy to demonstrate that if the atomic wavefunctions are 
the eigenstates of the parity and time-reversal operators, 
the induced magnetic moment vanishes. However, in the 
presence of the CP-odd interactions, V CP , there appears 
a tiny E- field-induced magnetic moment. To emphasize 
the essential role of CP-violation in the generation of the 
magnetic- moment, we will use CP superscript with the 
magnetic moment, (/i. CP ). The interaction V CF can be 
due to electron EDM or CP-odd weak neutral currents, 
and we will specify the particular forms of V CF in See- 



In these formulas, the summations are carried out 
over the eigenstates of the atomic Hamiltonian H ai 
Ha^p 1 " 1 ) = E p \^^). The derived third-order expression 
can be presented in a more compact and symmetrical 
form using the resolvent operator 1Z = (Eq — H a ) , 

(/x CP ) =2(0|U CP ft//ftV cxt |0)+ (7) 
2{o\ fJ ,nv CP nv cxt \o) + 2{0\ t inV c * t KV CP \0) . 

The three above contributions differ by permutations of 
the operators fi, V CP and V cxt . 



B. Dirac-Hartree-Fock approximation 

Having derived a general third-order expression for the 
induced magnetic moment, Eq. J7J), here we proceed with 
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the atomic-structure part of the evaluation. We em- 
ploy the conventional Hartree-Fock (HF) or independent- 
particle approximation for that purpose. In this ap- 
proach, the atomic many-body wavefunction is repre- 
sented by the Slater determinant composed of single- 
particle orbitals. These orbitals are determined from a 
set of the HF equations. Using a complete set of Slater 
determinants, the contributions to the induced magnetic 
moment, Eq.^BJ), may be expressed as 



<m cp > 1 =2£ 



V CP u V ext 

an H"nm * ma 
{£■171 &a n £a) 



(8) 
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(10) 

(11) 
(12) 
(13) 



Here indexes a and b run over single-particle orbitals oc- 
cupied in |^o), indexes m and n run over virtual orbitals, 
and Si are the energies of the HF orbitals. 

It is well known that the relativistic effects are essential 
for the non-vanishing contributions to energy levels due 
to EDMs (Schiff theorem). Moreover, in Section llll Al we 
will demonstrate that the relativism enters into the cal- 
culations of CP-violating polarizability in the enhanced 
fashion: one also needs to incorporate relativistic correc- 
tions to electric- and magnetic-dipole matrix elements 
and energies entering Eq. 10}© ■ We include the relativis- 
tic effects by directly solving Dirac-Hartree-Fock (DHF) 
equations 

(c(a • p) + /3c 2 + Kuc + Vdhf) u t {r) = e t u t {r) , (14) 

where Vn UC is a potential of the Coulomb interaction with 
a finite-size nucleus and Vdhf is non-local self-consistent 
DHF potential. 



C. Matrix elements 

We use the following ansatz for the Dirac bi-spinor 
1 ( iP nK [r) H rem (f) 



r I Qn K (r) fi_« m (r) 



(15) 



where P and Q are the large and small radial components 
respectively and Q is the spherical spinor. The angular 
quantum number k = (I — j) (2j + 1). 



In particular, the reduced matrix elements of the 
magnetic-dipole and electric-dipole moment operators 
between two bi-spinors are given by 



HHI&) 



(a\\D\\b) 



1 



(n a + K b ) (-K a \\Ci\\Kb) x 



(16) 



rdr{P a (r)Q b {r)+Q a (r)P b (r)}, 

n 

K«||Ci||K 6 ) X (17) 

X) 

r dr{P a (r)P b (r)+Q a (r)Q&(r)}, 



Ci(f) being normalized spherical harmonic. 

At this point we would like to specify particular forms 
for the CP-odd interaction V CP . We will distinguish 
between the electron EDM coupling T/ CP > EDM and weak 
neutral-current (NC) interactions 1/ CP > NC . The EDM in- 
teraction of an electron with an electric field Ew^ can be 
written in four-component Dirac notation as [13|: 



y-CP,EDM = 2 , 




er • £- mt 



(18) 



The matrix element of this interaction is given by 



V. 



CP, EDM 



d e ^2Z J ^Qa(r)Qb(r)}s Ka 



where we assumed that the dominant contribution is ac- 
cumulated close to the nucleus (of charge Z) so that £- lrA 
can be approximated by the nuclear field. The selection 
rules with respect to angular quantum numbers m and k 
arise because V CP is a pseudoscalar. 

Recently there was a renewed interest to CP-odd weak 
neutral current interactions of electrons with nucleons 
[l2| . It is known that in atomic experiments EDM of 
the electron is indistinguishable from the scalar CP-odd 
weak neutral currents 



LgUlS. 

m 



V 



CPNC 



. Gf 
.G F Z 



(Zk[ + ATfeD7075P(r) 
^r c 7075P(r), 



(20) 



where Gf = 2.2225 x 10 a.u. is the Fermi constant, 
fc^' n are dimensionless constants of the scalar P, T-odd 
weak neutral currents for proton and neutron (kf uc = 
+ ^fc"). Further, Z and N are the numbers of protons 
and neutrons in the nucleus, 70,5 are Dirac matrices, and 
p(r) is the nuclear density. 



III. RESULTS FOR RARE-GAS ATOMS 

The derived HF expressions hold for any atomic or 
molecular system with a state composed from a single 
Slater determinant. Below we will carry out calculations 
for the rare-gas atoms He through Rn. These closed-shell 



4 



atoms have a ^So ground state and, due to the spheri- 
cal symmetry, the CP-violating polarizability is a scalar 
quantity, i.e., the induced magnetic moment is parallel to 
the applied electric field. The intermediate many-body 
states in Eq. ar< 3 particle- hole excitations, with the 

total angular momenta of J — or J = 1 , depending on 
the multipolarity of the involved operator. 

To carry out the numerical evaluation, we solved the 
DHF equations in the cavity using a B-spline basis set 
technique by Johnson et al. 0]. The resulting set of 
basis functions is finite and can be considered as numer- 
ically complete. In a typical calculation we used a set 
of basis functions expanded over 100 B-splines. An ad- 
ditional peculiarity related to the Dirac equation is an 
appearance of negative energy states (e m < —m e c 2 ) in 
the summation over intermediate states in Eq. (|5|l- (|13[) . 
In our calculations we used the so-called length-form of 
the electric-dipole operator, Eq. I|17fl and we found the 
contribution of negative-energy-state to be insignificant. 



Atom 


Z 


/3 UP /d e 




He 


2 


3.8[-9] 


2.4[-22] 


Ne 


10 


2.2[-6] 


1.5 [-19] 


Ar 


18 


7.4[-5] 


5.2[-18] 


Kr 


36 


3.6[-3] 


3.1[-16] 


Xe 


54 


4.5[-2] 


5.3[-15] 


Rn 


86 


1.07 


2.2[-13] 



TABLE I: CP-violating polarizability, j3 , in Gaussian 
atomic units, for rare-gas atoms. CP-violation is either due 
to the electron EDM, d e , or due to the neutral currents 12011 . 
Notation x[y] stands for x x 10 v . 
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FIG. 1: (Color online) Dependence of the CP-violating polar- 
izability /3 CP on the nuclear charge Z for rare-gas atoms. CP- 
violation is due to the electron EDM, d e . The ratio /3 CP /d e 
is given in atomic units. 

Numerical results for rare-gas atoms are presented in 
Table H] and also plotted in Fig. ^ In Table [I] the values 



in the column marked /3 CP jd e were computed directly, 
while the values (3 CP /fc" uc (the last column) were ob- 
tained from (3 CP /d e as explained in Section lill Al 

From Fig. ^ we observe a pronounced dependence of 
the values on the nuclear charge Z. Such a steep scal- 
ing of the CP-odd polarizabilities is expected from the 
considerations presented in Sec. IIII Al 

To illustrate the (doubly) relativistic origin of the CP- 
odd polarizability (3 CF , we compile values of various con- 
tributions to /3 CP in TablellDfor an isolated Xe atom. Ap- 
parently, the dominant contributions are from (/i CP )i ja , 
Eq. ©, and (/i CP )i,f,, Eq. ©, but there is strong can- 
celation between these two terms. As we will see below, 
this cancelation is not accidental. 



k 






«/4 




sum 




1 


-0.108 




0.132 




2.44[- 


2] 


2 


6.53[- 


3] 


-6.63[- 


•5] 


6.46[- 


•3] 


3 


8.19[- 


3] 


5.13[- 


-3] 


1.33[- 


2] 


total 










4.42[- 


2] 



TABLE II: Contributions to CP-violating polarizability, 
/3 CP /d e , in Gaussian atomic units, for an isolated Xe atom. 
Each contribution is defined via Eq. ©-|H as ffi p a /d e = 
{/Jp P )k,a/(d e £o)- CP-violation is due to the electron EDM, 
d e . Notation x[y] stands for x x 10 y . 



A. Z 5 scaling and relation between EDM and NC 
contributions 

Let us consider non-relativistic limit of Eqs. (0]-[HJ). 
The magnetic moment operator is reduced to the form: 

M=f(2* + 1)- (21) 

This operator can not change electronic principal quan- 
tum numbers. Because of that the contributions © and 
© vanish, as there // should mix occupied and excited 
orbitals. Thus, we are left with the single term which 
can be further split in two parts (jHJ and 10 . We will show 
now that these two parts cancel each other. 

Indeed, in the non-relativistic approximation the oper- 
ator V is given by a scalar product of the spin vector 
and the orbital vector. Therefore, in the LS'-coupling 
scheme it can couple the ground state only with ex- 
cited states 3 Po- Operator (|21|l is diagonal in the quan- 
tum numbers L and S and can couple 3 Pq only with 3 Pi . 
To return back to the ground state, the dipole operator 
ycxt nag £ COU pi e 3p i w j^] 1 15^ However, this matrix 

element vanishes in the non-relativistic approximation. 
The above states 3 Pq } i are formed from the excited elec- 
tron and a whole in the core, which correspond to two 
expressions (|SJ) and @ . We conclude that these two con- 
tributions exactly cancel in the non-relativistic approxi- 
mation. 
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The matrix element ^PiIV^^Sq) is proportional to 
the spin-orbit mixing, which is of the order of (aZ) 2 . It 
follows from 1)16(1 that relativistic correction to operator 
((21(1 is of the same order. This correction accounts for the 
nondiagonal in the principle quantum numbers matrix 
elements of fi and leads to the nonzero values of the terms 
(J3J) and JHl- Thus, we see that all three terms in J3J are 
suppressed by the relativistic factor (aZ) 2 , in agreement 
with numerical results from Table ITT1 

Matrix elements of the CP-odd interaction V CP de- 
pend on the short distances and rapidly decrease with 
quantum number j. To a good approximation it is pos- 
sible to neglect all matrix elements for j > 3/2. For the 
remaining matrix elements between orbitals S1/2 and P1/2 
an analytical expression can be found in |l3[ : 

i a n 2 7 3 pEDM 
/„ |t/CP,EDMi \ lOO^ft . . 

{si/2\V \p 1 / 2 ) = - — 7 \TF?~ e > i 22 ) 



^l/ 2 |^ CP ' NC b 1/2 ) 



3 {vsVpY' 2 
G F aZ 3 R NC 



(23) 



where we use effective quantum numbers v = (— 2e) 1 / 2 . 
R EBM and i? NC are relativistic enhancement factors: 



I? 



EDM 



7(47 2 



fi Nc = «r^j;: 5 , Z=54 



1, Z = l, 

1.4, Z = 54, (Xe), 

2.7, Z = 86, (Rn), 

1, 



r 2 (2 7 + i) 



1.7, Z = 86, 



(24) 
(25) 



where 7 = \/l — (aZ) 2 and the radius of the nucleus is 
taken to be r N = 1.2 (Z + N)^ 3 im. 0. 

We see that both CP-odd operators scale as Z 3 R with 
relativistic enhancement factors R given by (|24|l and (1251) . 
This scaling adds up with relativistic suppression (aZ) 2 
discussed above to give overall scaling Z 5 R. This scaling 
agrees with our numerical calculations and Fig. ^ 

Because of the similarity between matrix elements 1(22(1 
and of operators y cp ' EDM and l/ Cp . NC , there is no 
need in calculating independently the NC contribution 
to (3 . It is sufficient to substitute matrix elements l|22f) 
in all equations with matrix elements ((23(1 . Comparing 
these expressions we find that to get the contribution to 
/3 CP induced by the CP-odd weak neutral currents we 
need to make following substitution: 



d e 
er 



dNC 

•^x^^edmC (26) 



where ro is the Bohr radius and i? JjjUM and R are given 
by (|231 and • The accuracy of Eq. is typically 
15 - 20%, which is sufficient for our purposes. It was 
used to calculate the last column of Table [I] 



IV. LIMITS ON ELECTRON EDM FROM 
MEASUREMENT OF CP-ODD 
POLARIZ ABILITY 

Here we envision the following experimental setup (see 
Fig. |2| to measure the CP-violating polarizability: A 
strong electric field £q is applied to a sample of diamag- 
netic atoms of number density n. A macroscopic mag- 
netization arises due to the CP-violating polarizability. 
This magnetization generates a very weak magnetic field 
B. One could measure this induced magnetic field and 
set the limits on the electron EDM or other CP-violating 
mechanisms. In particular, for a spherical cell the maxi- 
mum value of the generated magnetic field at the surface 
of the sphere can be related to the CP-violating polariz- 
ability as 



Bmax — — n (3 CP £ . 



(27) 



Clearly, one should increase the number density to en- 
hance the signal, and it is beneficial to work with a dense 
liquid or solid sample. 



Magnetometer 




FIG. 2: (Color online) A scheme for measuring CP-violating 
polarizability. 

Among the rare-gas atoms, considered here, xenon has 
the most suitable properties for such an experiment: Xe 
is the heaviest non-radioactive atom, it has a large num- 
ber density (n ~ 10 22 1/cm 3 ), and liquid Xe has a high 
electric field breakdown strength (£0 ~ 4 x 10 5 V/cm). 
Our calculations in Section IlIII were carried out for iso- 
lated atoms. However, in a liquid, there are certain en- 
vironmental effects (such as confinement of electronic 
density) that affect the CP-violating signal. To esti- 
mate the confinement effects in the liquid, we employ 
the liquid-cell model. The calculations are similar to 
those presented in Ref. In brief, we solve the DHF 
equations for a Xe atom in a spherical cavity of radius 

1/3 

Rcav = (37 ^7) > with certain boundary conditions im- 
posed at the cavity surface. For a density of LXe of 
500 amagat 01 1 ^cav — 4.9 bohr. For a solid state, 
i? cav ~ 4.2 bohr and we use the latter in the calcula- 
tions (see discussion in Ref. 0])- Technically, we applied 
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the variational Galerkin method on a set of 100 B-spline 
functions |15| . We find numerically that compared to 
an isolated atom, the CP-violating polarizability of a Xe 
atom in LXe is reduced by about 65%, 

/3 CP (LXe) w 1.5 x 10~ 2 4 . (28) 

From Eq. (|27|l it is clear that the more sensitive the 
measurement of the B-field, the tighter the constraints 
on (3 CP (and d e ) are. Presently, the most sensitive 
measurement of weak ma gnet ic fields has been carried 
out by Princeton group |18| . Using atomic magne- 
tomctry, this group has reached the sensitivity level of 
5.4 x 10~ 12 G/VHz. The projected theoretical limit [hJ 
of this method is 10 _13 G/\/Hz- Notice that this esti- 
mate has been carried out for a sample of volume 0.3 
cm 3 . According to Romalis 0], the sensitivity increases 
with volume V as U 1 / 3 , so a 100 cm 3 cell would have 
an even better sensitivity of about 10~ 14 G/Hz 1 / 2 . More 
optimistic estimate, based on nonlinear Faraday effect in 
atomic vapors [2(|, is given in Ref.[2l|; here the projected 
sensitivity is 3 x 10~ 15 G/a/Hz. 

Assuming 10 days of averaging, the most optimistic 
published estimate of the sensitivity to magnetic field [2ll ] 
leads to the weakest measurable field of B ~ 3 x 10 -18 G. 
Combining this estimate with the breakdown strength of 
the E-field for LXe, £q ~ 4x 10 5 V/cm, and our computed 
value of CP-odd polarizability, Eq. (|28|l , we arrive at the 
constraint on the electron EDM, 

de(LXe) < 6 x 10~ 26 e • cm. (29) 

This projected limit is more than an order of magnitude 
worse than the present limit on the electron EDM from 
the Tl experiment 0, 4 (Tl) < 1.6 x 10~ 27 e • cm. It 
is worth emphasizing that the above limit has been ob- 
tained using B-field sensitivity estimate from Ref. |2lj : 



with the present sensitivity record |18|, the constraints 
of electron EDM are several orders of magnitude weaker. 
In other words, we find that a substantial improvement 
in the experimental sensitivity to weak magnetic fields 
is required before the CP-violating polarizability of LXe 
can be used for EDM searches. 



V. CONCLUSION 

To summarize, we have computed novel CP-violating 
atomic polarizabilities 0, f3 CP , for rare-gas atoms. We 
have derived third-order expressions for (3 and em- 
ployed the Dirac-Hartree-Fock method to evaluate the 
resulting expressions. We have elucidated the doubly 
relativistic origin of the polarizability and demonstrated 
strong Z 5 dependence on the nuclear charge. Finally, we 
evaluated a feasibility of setting a limit on the electron 
EDM by measuring CP-violating magnetization of liquid 
Xe. We found that such an experiment could provide 
competitive bounds on electron EDM only if the present 
level of experimental sensitivity to ultra-weak magnetic 
fields is improved by several orders of magnitude. 
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